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Abstract 

Smoke aerosols from biomass burning are an important component of the global 
aerosol cycle. Analysis of Aerosol Robotic Network (AERONET) retrievals of size dis- 
tribution and refractive index reveals variety between biomass burning aerosols in 
5 different global source regions, in terms of aerosol particle size and single scatter 
albedo (SSA). Case studies of smoke transported to coastal/island AERONET sites 
also mostly lie within the range of variability at near-source sites. Two broad ‘fami- 
lies’ of aerosol properties are found, corresponding to sites dominated by boreal forest 
burning (larger, broader fine mode, with midvisible SSA ~0.95), and those influenced 
10 by grass, shrub, or crop burning with additional forest contributions (smaller, narrower 
particles with SSA ~0.88-0.9 in the midvisible). The strongest absorption is seen in 
southern African savannah at Mongu (Zambia), with average SSA ~0.85 in the midvis- 
ible. These can serve as candidate sets of aerosol microphysical/optical properties for 
use in satellite aerosol optical depth (AOD) retrieval algorithms. The models presently 
is adopted by these algorithms over ocean are often insufficiently absorbing to repre- 
sent these biomass burning aerosols. A corollary of this is an underestimate of AOD 
in smoke outflow regions, which has important consequences for applications of these 
satellite datasets. 


1 Introduction 

20 For several decades, satellite observations have provided a powerful tool for monitoring 
many aspects of the Earth system, including the atmospheric aerosol loading. Quan- 
tities such as aerosol optical depth (AOD) have generally been retrieved with lower 
uncertainties over oceans than land surfaces, due to the comparative homogeneity 
of open ocean surface properties and general lack of strong oceanic aerosol point 
25 sources. Despite this, significant differences can still exist between AOD retrieved us- 
ing different instruments or algorithms, in both clean and polluted conditions, and for 
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real measurements as well as simulated data (Remer et al., 2008; Kokhanovsky et al., 
2010; Sayer et al., 2012a). 

To some extent, these differences can be the result of differences between sensor ra- 
5 diometric calibration, cloud screening, or sampling-related issues such as pixel/retrieval 
weighting and averaging techniques (Mishchenko et al., 1999; Kaufman et al., 2005; 
Kahn et al., 2007; Levy et al., 2009; Sayer et al., 2010b). However, even over oceans, 
measurements made by past and current passive satellite sensors generally do not 
provide sufficient information to retrieve unambiguously all relevant parameters of in- 
io terest, such as spectral and directional surface reflectance, and aerosol microphysical 
properties and vertical distribution (Hasekamp and Landgraf, 2005). Thus, algorithms 
must make assumptions about these quantities, tied to the particular strengths and 
weaknesses of the instrument in question. Specifically, most algorithms parametrise 
oceanic surface reflectance as a combination of wind-roughened Sun-glint and foam, 
is with an additional contribution linked to oceanic chlorophyll concentration (e.g. Sayer 
et al., 2010a). Additionally, aerosols are often assumed to consist of a mixture of 
components whose total abundance and relative weight are varied in order to best 
match the top-of-atmosphere (TOA) radiance observed by the sensor (e.g. Martonchik 
et al., 1998; Mishchenko et al., 1999; Remer et al., 2009; Thomas et al., 2009; Sayer 
20 et al., 2012a). These algorithmic assumptions also contribute to uncertainty in retrieved 
aerosol properties, and the disparity between available satellite-derived datasets. 

As the AOD increases, so does its contribution to the TOA signal in the shortwave 
(and in some cases longwave) spectrum, and so the sensitivity of the retrieval algorithm 
to these assumptions about aerosol composition (e.g. Bulgin et al., 2011; Lee et al., 
25 2012). It is therefore important that these assumptions are realistic, in order to minimise 
errors in retrieved AOD. Existing datasets typically use microphysical properties derived 
from in-situ measurements, semi-empirical considerations, or else leverage retrievals 
of microphysical properties from sources with a higher information content, such as 
those in the Aerosol Robotic Network (AERONET; Holben et al., 1998). An excep- 
tion to this is the Polarization and Directionality of the Earth’s Reflectance (POLDER) 
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sensor, whose multidirectional and polarisation measurements offer an increased in- 
formation content, allowing algorithms to use weaker constraints about microphysical 
assumptions as compared to other sensors (Dubovik et al., 2011; Hasekamp et al., 

5 2011 ). 

Aerosol over the ocean is typically composed of a combination of hydrated sea 
salt particles and biogenic organic aerosols (O’Dowd and de Leeuw, 2007), with 
(regionally- and seasonally-dependent) contributions from transported continental 
aerosols (such as sulphates, nitrates, carbonaceous aerosols from industry or biomass 
io burning, and mineral dust). Sayer et al. (2012b) present an AERONET-derived model 
for unpolluted and dust-free marine aerosol. Microphysical properties of mineral dust 
(McConnell et al., 2008; Ansmann et al., 201 1 ; Su and Toon, 201 1 ; Lee et al., 201 2; Ry- 
der et al., 2013) and continental/industrial pollution (Dubovik et al., 2002; Aiken et al., 
2009; Giles et al., 2011, 2012; McMeeking et al., 2011) have been well-studied using 
is both AERONET data and field campaigns or monitoring sites. The focus of this study is 
to examine the microphysical properties of smoke-dominated aerosol mixtures trans- 
ported to oceanic regions, and compare these to the properties of smoke nearer to 
source regions. 

Biomass burning is an important part of the global aerosol burden. Smoke aerosols 
20 near their source regions are strongly optically-dominated by fine-mode absorbing par- 
ticles, with properties dependent on the substance which is burning as well as type 
of combustion (flaming vs. smouldering) and temperature/moisture content (e.g. re- 
gional/global reviews by Streets et al., 2003; Reid et al., 2005a, b; Janhall et al., 2010). 
Comparatively less well-studied are properties after medium- and long-range trans- 
25 port to the oceans, where processes such as chemical ageing, wet/dry deposition, and 
mixing with air masses of different origins may alter the microphysical properties of 
the columnar aerosol. Note, however, that some chemical changes within the aerosol 
occur fairly rapidly after emission (such as coagulation and moisture uptake), so even 
for near-source sites time differences of order of an hour or less can mean aerosol 
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microphysical and radiative properties differ from those at the time of emission (Radke 
et al., 1991 ; Hobbs et al., 1997; Martins et al., 1997; Magi and Hobbs, 2003). 

The analysis herein proceeds as follows. Section 2 discusses the AERONET data 
and assorted notation used. Then, the microphysical properties of smoke aerosols 
from AERONET sites near major global biomass burning source regions are examined 
in Sect. 3. As there are limited in situ measurements or AERONET sites located in the 
common biomass burning oceanic outflow regions, these climatological properties are 
then compared to cases of occasional smoke transport to coastal/island AERONET 
sites in Sect. 4. Section 5 illustrates potential biases in existing over-ocean satellite 
AOD datasets as a result of insufficient absorption in presently-assumed aerosol mi- 
crophysical properties, and Sect. 6 provides a perspective on the importance of the 
results. 


2 AERONET data 

The sun photometers used by AERONET measure spectral direct-beam solar radi- 
ance, as well as directional diffuse radiance in the solar almucantar. The former are 
used to determine columnar spectral AOD and water vapour, provided at a temporal 
resolution of approximately 10-1 5 min. Throughout this study, the AOD is denoted i x 
where A is the wavelength in nm. AERONET direct-Sun AOD has a typical uncertainty 
of 0.01-0.02 (Holben et al., 1998; Eck et al., 1999) and is provided at multiple wave- 
lengths (dependent on site) between 340 nm and 1640nm. The Angstrom exponent, 
a = -d Inr/c/lnA, is also derived by AERONET; references to a herein indicate that 
calculated over the wavelength range 440 nm-870 nm (by regression of all available 
over that wavelength range). 

The spectral AOD is also used, along with almucantar measurements of sky radiance 
taken over a large range of scattering angles, in the inversion algorithm of Dubovik and 
King (2000), Dubovik et al. (2006). This retrieves aerosol volume size distribution (in 
22 logarithmically-spaced bins with radii r from 0.05 pm-15 pm) and spectral complex 
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refractive index at 440 nm, 675 nm, 870 nm, and 1020nm. Cloud screening and other 
quality checks (Smirnov et al., 2000, Holben et al., 2006) are performed in order to 
remove potentially unreliable retrievals. Retrievals passing these checks are denoted 
“level 2”; for biomass burning aerosols and moderate aerosol loadings (t 440 ~ 0.4), as- 
sociated uncertainties for such retrievals are 25% on the binned size distribution (for 
0.1 pm</'<7pm, and larger on the tails), 0.04 on the real part of the refractive in- 
dex, and 30% on the imaginary part of the refractive index, giving an uncertainty in 
aerosol single scattering albedo (SSA, denoted ® 0 ) of approximately 0.03 (Dubovik 
et al., 2000). Higher AOD can decrease some of these uncertainties further. Size dis- 
tributions remain reliable at lower aerosol loadings, although refractive index uncertain- 
ties increase significantly. 

The notation adopted herein for aerosol microphysical properties follows that of 
Sayer et al. (2012b) and numerous other AERONET-based studies. The number size 
distribution d N(r)/d\u(r) describes the number of aerosol particles with radius in the 
infinitesimal size range r ± d ln(r); the related volume size distribution is 

dV(r) _ 47zt 3 d N{r) 
d\n{r)~ 3 d\n(r)’ 

for spherical particles. Total columnar aerosol particle number (C n ) and volume (C v ) 
are obtained by integrating these distributions over ln(r). The logarithmic volume mean 
radius (r v ) is a frequently-used metric of average aerosol particle size, defined 
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also often-used is the effective radius (r eff ), which is the ratio of the third to the second 
moment of the size distribution. The broadness of the distribution is often characterised 
by its spread (also called width), o, where 



The geometric standard deviation, e a in this notation, is sometimes used instead of 
o. The above definitions are independent of the shape of the size distribution. In prac- 
tice, aerosol size distributions are often represented as a combination of lognormally- 
distributed components (because real size distributions often closely follow a lognor- 
mal distribution, as well as for the mathematical/computational conveniences afforded 
by these distributions). In this case, the number size distribution is defined as a sum- 
mation over n c components by 


1 /\n{r)-\n(r nJ ) 
dm = ^ _^ 0 -2 { — o, — 
dW) £ j2Hoi 


(4) 


and the modal radius for each component is also its median and geometric mean. The 
equivalent formulation for aerosol volume is arrived at by substituting r n with r v , and 
C n with C y . For a lognormal component, conversions between number and volume 
quantities are (Sayer et al., 2012b) 

-,3(J 2 


r v = r n e 


C v = ? r n 3 e 4 - 5 ° 2 C n 


(5) 

( 6 ) 
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and additionally 
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Observed aerosol size distributions are typically (although not exclusively) bimodal 
(n c = 2), e.g. Dubovik et al. (2002). In these cases, and adopted here, the smaller (fine) 
mode properties are denoted with a subscripted f (i.e. r vf , o f ) and the larger (coarse) 
mode properties with a subscripted c. Analagously, fine and coarse mode AOD are de- 
noted r f A and r c A respectively. This AERONET dataset defines fine and coarse mode 
properties by locating the inflection point in the retrieved size distribution; as the distri- 
butions studied in this work consist of two distinct peaks, results are numerically only 
weakly sensitive to the precise fine/coarse demarcation point (coarse mode more so 
than fine, when r c is low). 

Propagating the aforementioned errors on retrieved aerosol size distribution 
(Dubovik et al., 2000) onto these distribution parameters (for a typical biomass burning 
model from Dubovik et al., 2002, and assuming errors in each bin of d V {r) / d \r\{r) are 
uncorrelated) suggest uncertainties on AERONET-derived r vf of 0.01 pm, o { of 0.06, r v c 
of 0.2 pm, and o c of 0.06 for individual retrievals. Throughout this work, “fine” refers to 
accumulation-mode aerosols; smaller “nucleation-mode” aerosols may also be present, 
with some of these existing on the smaller-/" tail of the fine mode, and others below this 
limit (although in that case so small as to be optically inactive in this spectral range). 
Mie theory is used throughout; although fresh smoke particles may be nonspherical, 
they rapidly coagulate to form aggregates whose radiative properties (in contrast to 
nonspherical coarse aerosols such as mineral dust) have been found to be modelled 
adequately by Mie theory (Kahn et al., 1 997; Mishchenko et al., 1 997; Reid and Hobbs, 
1998; Dubovik et al., 2000, 2006; Reid et al., 2005b). 
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3 Properties of smoke aerosol from near-source sites 
3.1 Data selection 

Level 2 version 2 AERONET inversion data (Holben et al., 2006) from eight AERONET 
sites, shown in Fig. 1 (and with key information and references in Table 1), are used 
to investigate the microphysical properties of smoke aerosols over a range of different 
burning regimes. Some of these sites were previously used in a similar analysis by 
Dubovik et al. (2002); the analysis here benefits from a decade more observations, 
as well as improvements made to the version 2 AERONET database since that time 
(e.g. surface reflectance inputs; Eck et al., 2008). Giles et al. (2012) also examined a 
few of these sites, although that study was focussed on characterising global aerosol 
absorption properties from a variety of types, rather than creation and examination of 
microphysical models to represent these aerosols. State/country information for each 
site is given parenthetically the first time a given site is discussed in the text. 

Some regions contain several sites in relatively close proximity; in this case, typically 
the site with the largest data record and smallest influence from other aerosol sources 
was chosen. Manual inspection reveals that data from these nearby sites are generally 
similar. For example, Mukdahan (Thailand) was chosen as a key site for south-eastern 
Asia, although quantitatively similar results are obtained if data from Pimai (also in 
Thailand) are used instead. Both sites primarily sample nearby agricultural burning, 
and forest burning from elsewhere in the region. In contrast, the site in the city of Chi- 
ang Mai (north-western Thailand) was not used as, while close to biomass burning 
source regions, it also has a significant urban contribution to the aerosol, coupled with 
the local topography leading to a “trapping” of pollutants (Gautam et al., 2012; Janjai 
et al., 2012). Another point to note is the omission of sites from the Sahel; although 
AERONET contains several long-term sites in this region, the dry season biomass 
burning peak here (November-February) due to burning of grassland, cropland, and 
shrubs coincides with strong dust activity, such that the aerosol is normally a mixture 
of biomass burning smoke and dust (Pandithurai et al., 2001; Roberts et al., 2009). 

25021 
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Johnson et al. (2008) and Eck et al. (2010) found cases where these smoke aerosols 
were highly absorbing, with midvisible SSA as low as 0.76, although the more com- 
5 mon cases of mixed smoke and dust in this region are less absorbing. Limiting to 
smoke-dominant cases leaves a small data volume (and an unrepresentative sample 
of the true nature of the aerosol in that region). Additionally, the smoke and dust in 
this area often occur in distinct vertical layers and so a microphysical model based 
on column-average aerosol properties, when applied to real satellite data, could po- 
io tentially be misleading or inaccurate (Kim et al., 2009). A similar rationale applies to 
sites in the Indo-Gangetic Plain and Himalayas (e.g. Pokhara, Nepal), as springtime 
biomass burning overlaps with transported dust, and some sites additionally have a 
significant urban component to the aerosol (Gautam et al., 201 1 ). 

At each site, in order to remove likely non-smoke cases from the data (as burning is 
is seasonal and aerosol of other types may be transported periodically to sites), two sets 
of filters are applied. Firstly, it is required that r 440 > 0.4 and a > 1, to ensure reliabil- 
ity of retrieved refractive index/SSA and remove potential cases dominated by mineral 
dust aerosol (which has smaller a, e.g. Eck et al., 1999). Secondly, data are restricted 
to the main burning season for those sites where this is well-defined (see Table 1); 
20 some sites, such as Mongu (Zambia), exhibit fairly constant burning during through- 
out a period of several months, while others, such as Bonanza Creek (Alaska, USA), 
have a comparatively low baseline AOD punctuated by episodes of burning (which 
tend to occur in local summer months in dry years). As biomass burning is the domi- 
nant reason for high-AOD conditions at all these sites, the main effect of these filters is 
25 to remove some outlying results. Additionally, data from the year 2008 were excluded 
from CUIABA-MIRANDA (Mato Grosso, Brazil; hereafter Cuiaba), and from 2002 and 
2004 from Jabiru (Australia), due to suspected instrumental issues which may affect 
retrieved SSA during these periods. 
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3.2 Overall and comparative microphysical/optical properties 
3.2.1 Overview of summary Figures and Tables 

The median retrieved volume size distributions corresponding to these biomass burning 
5 cases at each site are shown in Fig. 2. These reveal the presence of a strong fine mode 
with peak radius in the range 0.1-0. 2 pm, and a secondary coarse mode with peak 
radius in the range 3-5 pm. Both modes appear approximately lognormal, although the 
coarse mode has a slight low-radius skew; the minimum in the overall aerosol volume 
size distribution is most commonly from 0.5-0. 9 pm. On an individual-retrieval basis, 
io and dependent upon the total AOD, the fine mode typically accounts for 40-90 % of the 
total aerosol volume (more at higher AOD) and 80-95% of the AOD at 550 nm. 

The AOD-dependence of aerosol size distribution parameters has previously been 
noted for a wide variety of aerosol types (e.g. Dubovik et al., 2002). Therefore, these 
data were also examined for these tendencies. An example of this process is shown in 
is Fig. 3 for fine-mode radius and spread (r vf , <j f ) at Alta Floresta (Mato Grosso, Brazil). 
To mitigate the effect of outliers, size distribution parameters more than two standard 
deviations away from the median at a given site were discarded for this portion of the 
analysis, and then a linear least-squares fit of parameter against fine-mode AOD per- 
formed. The exclusion removed typically ~5% of data, and did not lead to a significant 
20 difference in regression statistics, at most sites. Resulting linear regression relation- 
ships are shown in Table 2 for r f 550 (and Table 3 for r f 440 , for reference), and illustrated 
later. The regression fits were performed against fine-mode AOD as this is likely a more 
reasonable metric of the contribution of biomass burning to the total aerosol burden. 
This is consistent with the primary emissions from biomass burning being small aerosol 
25 particles and aerosol precursor gases. For all sites the coarse-mode AOD was typically 
~0. 02-0. 04 and roughly independent of the total AOD. Flence, similar results are ob- 
tained if total AOD is used instead. For a similar reason (low coarse-mode AOD and 
noisier size distribution parameters), no regression was performed for coarse-mode 
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size distribution parameters as a function of coarse AOD; instead, median values of 
r v c and o c are given in Tables 2 and 3. 

Aerosol refractive index was also found to be largely independent of AOD. The site- 
median refractive index, and SSA calculated using the relationships in Table 2 for the 
case of r f 550 = 0.5, t c 550 = 0.03 at each site, are given in Table 4. Although the AOD- 
5 dependence of model parameters means that SSA will itself be a factor of fine-mode 
AOD, the AOD-dependence of SSA is small according to these relationships (generally 
becoming less absorbing by ~0.01 as fine-mode AOD increases from 0.2-2). 

The relationships between all these properties and total columnar water vapour were 
also explored, but in most cases omitted for brevity, due to a lack of observed covari- 
io ability. This is consistent with the idea that biomass burning aerosols show weaker 
hygroscopicity than other fine-dominated aerosols (e.g. Reid et al., 2005a), and that 
ageing processes may have occurred rapidly and prior to the aerosols reaching the 
AERONET site. This could also be linked to the small gradient of the observed r f 550 - 
to-r vf relationships at these sites, which are up to about an order of magnitude weaker 
is than observed for more hygroscopic urban/industrial aerosol particles (e.g. Dubovik 
et al., 2002). Despite this, Reid et al. (2005a) did note that some measurements of 
biomass burning aerosols in scattered regions showed larger hygroscopic growth fac- 
tors than expected, although the reasons for this were not understood. Additionally, as 
the AERONET aerosol and water vapour data represent column-averaged quantities 
20 they do not provide information about the extent of vertical overlap between aerosol 
and moisture. Changes due to moisture uptake may also be masked amongst variabil- 
ity of aerosol properties from other sources. One corollary of this is that, when using 
these parameters as a source of biomass burning aerosol microphysical properties for 
satellite AOD retrieval, it may not be necessary to account for hygroscopic growth of 
25 these particles for some burning types. 

To facilitate a comparison between the different sites, Fig. 4 compares calculated 
size distributions and optical properties for the aforementioned case of r f 550 = 0.5, 
t c,55o = 0-03. Also shown are properties for the ‘fine-dominated’ aerosol model used 
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in the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Ocean Aerosol Retrieval 
(SOAR, Sayer et al., 2012a) algorithm, which is similar (/- eff f within 0.03 pm, SSA lower 
by ~0. 02-0. 03) to models used in the operational over-ocean Moderate Resolution 
Imaging Spectroradiometer (MODIS; Remer et al., 2009) algorithm, and within the 
range of aerosol components included in the Multiangle Imaging Spectroradiometer 
(MISR; Kahn et al., 2010) algorithm. 

3.2.2 Discussion of aerosol properties 

Considering the above data, the eight sites span a range of size distributions and 
optical properties but there are some similarities between them. Bonanza Creek and 
Yakutsk (Siberia, Russia), both boreal forest/peat burning sites (likely mostly dominated 
by smouldering combustion), show the weakest and roughly spectrally-neutral absorp- 
tion (® 0 ~0.95), and are among the sites with a larger volume radius and broader 
distribution. Linear relationships between aerosol properties and AOD or water vapour 
at this site show strong scatter, but if plotted together as a function of day of year 
(Fig. 5), some seasonal tendencies become visible. Both r vf and a f tend to be larger 
for higher r f 550 , which occurs most often in August at Bonanza Creek and more spo- 
radically at Yakutsk. This is likely to soil drying through the season making peat soils 
more flammable through time (Eck et al., 2009, and references therein). There is a 
suggestion of a bowl-shaped profile of r vf , with the lowest values from May-July and 
higher values earlier and later in the season. Columnar water vapour begins to increase 
around the start of June, and falls off during September. The SSA of the smoke is gen- 
erally higher from mid-May onwards. It is possible that these seasonal changes reflect 
changes in moisture and/or vegetation phenology, resulting in a change of the nature 
of burning. Strong scatter is expected at these sites in particular as burning can occur 
across wide tracts of the land at these latitudes and be transported long distances, i.e. 
smoke sampled at these two sites may be of diverse origin and age. 

In contrast, Cuiaba, Skukuza (South Africa), and Jabiru (grass, crop, and shrublands) 
are very similar to each other, with narrower lower-radius distributions, and are more 
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strongly absorbing (® 0 ~0.88-0.9 in the visible, and ~0.85 in the near-infrared). Note 
that Cuiaba can also sample forest burning from the north, while Skukuza includes 
a sulphate contribution from industrial sources (Piketh et al., 1999). Mukdahan has a 
similar SSA to these sites, although size distribution shape and asymmetry parameter 
closer to Yakutsk. For this set of sites relationships between aerosol properties and 
5 water vapour were not evident, except at Mukdahan, where a weak increase (~0.03) in 
SSA from February to April was mirrored by an increase in water vapour during these 
inversions from ~2.5-4cm. However this is likely linked to transport of air masses in- 
cluding pollutants from China/India, and an increase of forest burning relative to agri- 
cultural burning during this period, rather than hygroscopicity. 

10 Intermediate between the properties of these two groups is Alta Floresta, which sam- 
ples tropical woody burning from the nearby area. Relationships between aerosol fine- 
mode size/spread and water vapour (not shown) were found to be weak. The SSA 
at Cuiaba was found to be ~0.05 higher in October than August, while at the same 
time monthly median water vapour increased from ~2-4cm; a similar change in water 
is vapour at Alta Floresta through this period was not associated with any SSA change at 
this site. Together with the lack of change of size distribution parameters, this is consis- 
tent with the change in SSA at Cuiaba being driven with an increase in air masses con- 
taining transported smoke from forested regions through the season, rather than hygro- 
scopicity. The overall lack of relationship with water vapour content is consistent with 
20 Schafer et al. (2008), who noted almost no dependence of AERONET aerosol proper- 
ties on water vapour in Amazonia, except for the very highest AOD cases (r 550 > 1 .5). 
It is also worth noting that the area around the Alta Floresta site has become less 
forested due to agricultural conversion since the site was originally set up; analysis of 
the time series of retrieved aerosols properties at this site (not shown) reveals some 
25 interannual variability but no apparent trends or secular changes. 

Aerosols at Mongu and Skukuza have similar size distribution properties, although 
the former is significantly more absorbing (SSA lower by about 0.05), likely due to 
the contribution from non-absorbing sulphates to the aerosol at Skukuza (Piketh et al., 
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1999). Eck et al. (2013) found an increase in SSA of about 0.1 through the burning sea- 
son at both of these sites (as well as inferred from AERONET at Etosha Pan, Namibia, 
and satellite observations), attributed to a likely decrease of black carbon content from 
May through to November due to sampling of fires from different fuel types. This tem- 
poral variability was also seen to manifest as temporal variability in bias of satellite- 
retrieved AOD. At both of these sites, water vapour was low (typically ~1-2cm) and 
showed no links with aerosol properties. 

The asymmetry parameter ( g ) is 0.68 ±0.02 at 440 nm for all models (Fig. 4), but 
diverges to the range 0.43-0.57 at 1020nm (due to the differences in fine/coarse AOD 
partition at longer wavelengths and fine mode aerosol size), again roughly along the 
lines of the groupings given previously. Similar values and spectral dependence for 
fine-dominated aerosols were found in previous analyses (e.g. Dubovik et al., 2002), 
while g for dust aerosols tends to show less spectral variability. In terms of spectral 
dependence of AOD, most models give about 1 40 % AOD at 440 nm relative to 550 nm, 
and 30% at 1020nm relative to 550 nm, yielding values of a across the midvisible in 
the range 1.8-2. The one outlier here is Bonanza Creek, for which this fine/coarse 
mode combination gives (7=1.43, due to the larger and broader fine mode relative 
to the other sites. Note that although this linear a formulation is useful, it is only an 
approximation and caution should be taken when either using a to extrapolate AOD, 
or when comparing a determined across different wavelength ranges (e.g. Eck et al., 
1999). 

Examining the regression relationships in more detail, the fine-mode volume median 
radius for a moderate T f550 ~0.3 is typically ~0. 14-0. 15pm for the aforementioned 
grass/shrubland sites and larger (0.1 5-0.2 pm) for the wood-burning sites. The gradient 
of the r f i550 -to-r Vif relationship is in the range 0.01 3-0.02 at most sites. Low to moderate 
correlations for these relationships (0.23 at Skukuza, up to 0.62 at Cuiaba) reflect both 
the low range of AOD spanned at some sites, possible contributions from other aerosol 
sources, and the fact that the scatter around these relationships (shown later) of order 
0.01-0.02 is a similar size to the estimated uncertainty on AERONET retrievals of r vf 
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(mentioned previously). Similarly gradients and correlations with r f 440 instead of r f 550 
are slightly lower as the same data and scatter are essentially being stretched across 
the AOD axis. Positive correlations between r f 550 and o f are also observed, although 
are weak, likely for the same reasons of large scatter and limited dynamic range. All 
the above relationships are statistically significant at the 90% level, except for r f vs. 
o f at Jabiru (although again, many relationships are numerically small). It is important 
to note that, although these relationships have been derived in a climatological sense, 
the scatter about them comprises a significant fraction of the observed variablility, due 
to both retrieval noise and sources of true natural variability. 

Figure 6 shows that there is a fairly strong positive relationship between retrieved r vf 
and Of (r = 0.75). This is true both across the ensemble of sites as a whole, as well 
as if sites are examined individually (not shown). As the uncertainty on each quantity 
can be significant compared to the data range, a bivariate linear fit (e.g. York, 1966) 
was used, with uncertainty on r vf of 0.01 and on o f of 0.06. If the data are expressed 
in terms of number mean radius rather than volume (i.e. r nf rather than r vf ; Eq. 5), an 
opposing tendency is revealed: larger number mean radii are associated with narrower 
distributions, although the correlation is weaker (-0.46). A weaker correlation could in 
part be due to a larger uncertainty on the tails of the retrieved aerosol size distribution 
(Dubovik et al., 2000), to which r nf is more sensitive than r vf . This anticorrelation was 
also observed across 36 observations of fresh and aged smoke collated by Janhall 
et al. (2010). This could be caused by an initial broad distributions of emitted small 
particles coagulating into a narrower distribution of larger particles. 


4 Smoke transport to coastal/island sites 

Occasional cases of transported smoke have been observed at a wide range of coastal 
or island AERONET sites. Figure 7 shows smoke from several wildfires in California 
(USA), which are not uncommon in Northern Hemisphere summer and autumn, blown 
into the Pacific Ocean. Although not a rare occurrence, these sites are nonetheless 
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included in the “transported” rather than “near source” category in this analysis as 
there are insufficient cases passing over the AERONET sites to create a meaningful 
climatology of properties. 

Figure 8 shows an example of more long-range transport, namely Amazonian smoke 
transported south through South America, eventually passing over the AERONET site 
5 in Buenos Aires (Argentina) and out into the southern Atlantic Ocean. This pathway 
is confirmed by the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT; 
Draxler and Hess, 1998) model back-trajectory, for air arriving at Buenos Aires for ~2 
km altitude (HYPSPLIT estimates of the altitude of this air mass near the source were 
3-4.5 km). Aerosol properties at the sites in this figure are shown in more detail in 
io Fig. 9: r 440 at Buenos Aires is low for most of the period, while at the other sites it 
is larger, approaching 4 at Alta Floresta, and more variable (Alta Floresta and Cuiaba 
lack observations prior to 18 August). High r 440 and a at these sites is consistent with 
biomass burning; the fine mode effective radius and SSA at Alta Floresta and Ji Parana 
(Rondonia, Brazil) track each other reasonably closely, while the smoke at Cuiaba is 
is slightly smaller and more absorbing. During 21-23 August, the plume is transported 
from the northern region to Buenos Aires (Fig. 8), and the SSA matches that at Alta 
Floresta/Ji Parana (~0.9 at 440 nm), from where the plume is observed to travel. The 
absence of SSA retrievals outside this period at Buenos Aires is due to the low AOD 
(Dubovik et al., 2000). Note that the more-absorbing Cuiaba site lies outside of the 
20 main path of this plume. The fine mode effective radius at Buenos Aires is slightly larger 
during this period than before or after, and also larger than the near-source sites, which 
may be coincidental or may result from mixing with another air mass during transport. 
Following this event, AOD at Buenos Aires returns to typical low levels. Towards the end 
of the month satellite images (not shown) reveal that smoke from the forested region is 
25 blown over Cuiaba; the SSA and fine mode radius at Cuiaba increase to more closely 
match Alta Floresta/Ji Parana. 

Unfortunately, the number of cases where a smoke plume is observed at an 
AERONET site near-source and is conveniently transported past multiple other sites 
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over the course of days, thus providing the ability to track the evolution of the plume 
with AERONET, is very limited. Thus the comparison herein is more of a categorical 
nature: comparing how the body of case studies at island/coastal sites compares with 
the climatology and variability of land sites. Table 5 presents the list of AERONET inver- 
sions corresponding to such transported smoke identified at coastal/island sites (Fig. 1 ) 
5 used in this study. These data were individually identified as likely having a signifi- 
cant contribution from transported smoke by examination of satellite images, HYSPLIT 
back-trajectories, news/government agency reports, and guided by previous studies 
(Queface et al., 2003; O’Neill et al., 2005; Eck et al., 2009; Qin and Mitchell, 2009; 
Castro Videla et al., 2013). 

io Retrieved aerosol size distribution parameters and SSA for these cases are shown 
in Figs. 10, 11, 12, and 13, alongside those of the near-source AERONET sites which 
provide the closest match to these inversions. Properties for the “pure marine” aerosol 
model of Sayer et al. (2012b) are also shown for reference. This ‘best match’ was 
allocated subjectively, by identifying the source site for which the coastal/island site 
is data in question lay most neatly within the cloud of points. For example, some cases at 
Buenos Aires are more similar to Alta Floresta, while others are more similar to Cuiaba, 
reflecting the different potential origins of the smoke. The main point is to illustrate 
that there is overlap between near-source and transported smoke aerosol properties, 
rather than to tie each individual case down to a specific burning type. Additionally, the 
20 fine mode radius/spread comparison is shown for the aggregate of these data points 
in Fig. 6. Similar relationships as for the near-source cases are observed, albeit with 
slightly different regression coefficients and weaker correlation, perhaps due to the 
more limited data range. 

This comparison reveals that the fine-mode aerosol properties for these cases of 
25 smoke transported to ocean/coastal sites, in almost all cases, fall within the range of 
variability of smoke properties for the near-source sites. Some additional variability is 
observed, which could be attributed to additional ageing of the smoke aerosol particles 
and/or combination with aerosol particles from different air masses. In particular, most 
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of the coastal sites are in or near urban areas for which some additional urban/industrial 
aerosol contribution would be expected, as well as a potential marine aerosol contribu- 
tion. Some of the matchups are not particularly surprising due to relative geographical 
proximity, e.g. Skukuza (South Africa) and Inhaca (Mozambique), or Jabiru and Darwin 
(Australia). Others, such as the similarity between properties at Cuiaba and the cases 
5 of Californian wildfires, are more noteworthy. The main conclusion from this is that mi- 
crophysical models based on fine-mode smoke aerosol properties near source regions 
are also representative of smoke transported over the ocean, and so suitable for use 
in satellite AOD retrievals over ocean. Coarse-mode properties show more diversity 
between sites (and also often differ from pure marine aerosol properties), although as 
10 the coarse mode contribution to total AOD is minor this is not likely to be a significant 
source of error. 

Despite the similarities identified between near-source and coastal/island trans- 
ported smoke, there remain some notable differences. Inversions at the COVE site 
(off the coast of Virginia, USA) were for a case of smoke from Canadian boreal for- 
15 est fires in July 2002 (O’Neill et al., 2005). In terms of fine mode radius and aerosol 
SSA, these show a similarity to the climatological properties of boreal forest burning 
at Bonanza Creek (and also Yakutsk; Figs. 10, 12). However, <j f was smaller for these 
cases at COVE than commonly-observed at either of these sites. Additionally, there 
is comparatively large scatter in the observed fine-mode microphysical properties at 
20 these sites. This suggests that there may be significant variability in aerosol properties, 
perhaps influenced by variability in local conditions at the time of burning. A similarly 
large scatter is seen for Bach Long Vy (Fig. 12), in the Gulf of Tonkin (off the coast of 
Vietnam). Here it is likely that there is some urban/industrial component to the aerosol 
from nearby cities in Vietnam, flowing out of the Red River delta; this could explain the 
25 higher values of r vf and a f at this site. Additionally, this site’s location makes it more 
susceptible to instrumental degradation (e.g. moisture deposition and corrosion due to 
sea salt) than other sites, meaning retrievals may be less reliable here. A similar urban 
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contribution could explain the larger r vf at Sevastopol (Ukraine) and Singapore than 
the comparison site of Yakutsk. 

The largest outlier is Ascension Island, in the tropical Atlantic, which samples air 
masses including mixed Saharan/Sahelian dust and smoke aerosols from November- 
February, and central African smoke from June-October (only this second period con- 
tributed to the cases shown here). At this site, a f and aerosol SSA are a very close 
5 match to data from Mongu, while coarse-mode properties are similar to the pure mar- 
itime case (Fig. 11). Flowever, r vf is about 0.02 pm larger at Ascension Island than 
Mongu. This set of aerosol parameters is not a close match to any of the other near- 
source sites studied. 

The Ascension Island site is fairly remote and in a harsh environment (exposed to 
io salt from breaking waves, which may deposit on the instrument), which can lead to 
instrument problems more frequently than at some other AERONET sites. The most 
common symptom of these problems is an anomalously low SSA (typically up to 0.1 
lower than expected). Errors in retrieved size distribution parameters (and direct-Sun 
spectral AOD) as a result of these problems are expected to be minor. For this rea- 
ls son, AERONET inversions from Ascension Island were subject to additional calibra- 
tion/contamination checks beyond those applied for AERONET Level 2.0 processing 
as an extra precaution, and the data shown are for only those time periods where there 
was no indication of problems. Therefore it is likely that the higher r vf at Ascension 
Island than Mongu is a real characteristic of the aerosol transported to this area, rather 
20 than an artefact. Potential reasons include, again, additional aerosol ageing in the time 
from emission to arrival at Ascension Island, as well as the possibility that properties of 
the freshly-emitted aerosol are different. Some air masses reaching Ascension Island 
pass over parts of Africa north of Mongu, which are more heavily forested (Roberts 
et al., 2009). Unfortunately, there are no AERONET sites in this region. 
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25 5 Implications for satellite AOD retrievals 

Figure 4 revealed that the “fine-dominated” aerosol microphysical model used in over- 
ocean SeaWiFS processing (Sayer et al., 2012a) is intermediate in strength of absorp- 
tion between the more weakly-absorbing boreal forest sites, and the more strongly- 
absorbing tropical forest and grass/shrubland sites. Therefore, use of this model to 
retrieve AOD from satellite measurements in these latter cases is likely to result in an 
underestimate of AOD (or conversely overestimate AOD for cases of weakly-absorbing 
5 smoke). The most absorbing fine-mode aerosol component in the widely-used MODIS 
operational processing over ocean is slightly less absorbing than this SeaWiFS model 
(Remer et al., 2009). This suggests that the over-ocean AOD from some regional 
smoke aerosols in these datasets, which are large and seasonally-repeating features 
in some parts of the world, but have been validated only sparsely due to a lack of 
io ground truth data, may be understimated. In contrast the MISR aerosol mixtures in- 
clude aerosols with a midvisible SSA down to about 0.8 (Kahn et al., 2010), which may 
be sufficient, if the retrieval algorithm is successful in choosing an appropriate mixture 
(although note again rare cases of even lower SSA; Johnson et al., 2008; Eck et al., 
2010). All algorithms also include effectively nonabsorbing aerosol models. Although 
is this discussion focusses on retrievals over ocean, it is worth mentioning for complete- 
ness that incorrect SSA also leads to biased AOD retrievals over land (e.g. Ichoku 
et al., 2003; Eck et al., 2013, for smoke examples). 

Turning to observations, Fig. 14 compares midvisible AOD from these satellite prod- 
ucts against AERONET data at Ascension Island. Data from November-February are 
20 excluded to minimise the contribution of transported dust, which happens in this sea- 
son. The satellite:AERONET matchup protocol is as in Sayer et al. (2012a), namely 
AERONET data are spectrally interpolated to 550 nm and averaged within ±30 min of 
the satellite overpass, and satellite data are averaged within ±25 km of the AERONET 
site and restricted to only those retrievals meeting the dataset creators’ recommended 
25 quality assurance flags. For t 550 > 0.3 points, AERONET-derived a was most com- 
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monly ~1-1.6, consistent with mixed marine and smoke aerosols. The positive bias of 
MODIS and MISR data in low-AOD conditions (0-0.2) was noted in previous studies 
(Kahn et al., 2010, Sayer et al., 2012c). However, all datasets exhibit a low bias in AOD 
in conditions of elevated AOD, most notably SeaWiFS (although SeaWiFS has negligi- 
ble bias in clean conditions while the others have a positive bias; the change in bias as 
AOD increases is more similar between all datasets). 

To test the effect of aerosol absorption on satellite measurements, the 6S radia- 
tive transfer code (Vermote et al., 1997) was used to simulate TOA reflectance at 
5 wavelengths used for AOD retrieval by these sensors, for a variety of geometries, 
over an ocean surface with 6 ms -1 wind speed. Three aerosol types were considered: 
strongly-absorbing aerosol using the model for Mongu (Table 2, ® 0 ~0.85 at 440 nm); 
moderately-absorbing aerosol using the fine-dominated model of Sayer et al. (2012a) 
(® 0 ~0.95 at 440 nm); and the pure marine model of Sayer et al. (2012b) (® 0 ~0.99 
io at 440 nm). Then, for each wavelength simulated, these latter two models were used 
to retrieve AOD (reported relative to 550 nm) in each band, taking the Mongu case as 
‘truth’, and so calculate the AOD retrieval error. Although this does not mirror how the 
individual satellite algorithms mentioned previously function, it does provide a direct 
comparative baseline of the sensitivity of each wavelength to the assumed strength of 
is aerosol absorption. 

Figure 15 illustrates the results of this test for a solar zenith angle of 45°, viewing 
zenith angle of 10°, and relative azimuth angle of 135°. Similar patterns are observed 
at other common Sun/sensor viewing geometries (not shown). This figure shows the in- 
crease of TOA reflectance with AOD, with the increase being less pronounced for more 
20 strongly-absorbing aerosols, leading to a low bias in retrieved AOD if the real aerosol is 
less absorbing than assumed. The difference is larger for shorter wavelengths, linked to 
the larger aerosol signal and increased Rayleigh-aerosol interactions. This lends sup- 
port to the interpretation of Fig. 14, and suggests that future versions of satellite AOD 
retrieval algorithms should include an analogue for these strongly-absorbing aerosol 
25 particles. 
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Interestingly, although the current MISR algorithm includes strongly-absorbing 
aerosol mixtures, the change of bias from low-AOD to high-AOD conditions is similar 
to that in other sensors. This suggests that these mixtures may not always be chosen 
when needed (and the selection of an appropriate aerosol microphysical model to use 
in a system with a limited information content is a difficult problem in itself). The biases 
in high-AOD conditions in Fig. 14, typically from -0.1 to -0.3 for an AERONET AOD of 
0.5, are similar in magnitude to those suggested in Fig. 15 as being due to insufficient 
absorption in the assumed aerosol properties. These results provide evidence that the 
assumed SSA may be to blame - although other aforementioned factors (e.g. surface 
reflectance, calibration, pixel selection) may also contribute. 


6 Perspective 

Biomass burning is one of the major contributors to the global aerosol burden, with 
both natural and anthropogenic sources. Analysis of AERONET retrievals of size dis- 
tribution and refractive index revealed considerable variety between biomass burning 
aerosols in different global source regions, both in the size and strength of absorption 
of fine-mode particles. An advantage of AERONET is the long-term nature of obser- 
vations at some sites, and consistency in observation, retrieval, and quality-assurance 
procedures between sites, making it a useful tool for comparisons of this type. Derived 
aerosol properties were found to be within a similar range to those observed in other 
studies by a variety of techniques. Additionally, case studies of transported smoke from 
a larger range of AERONET sites also frequently fell within the range of variability for 
these near-source sites. 

Two broad “families” of aerosol properties were found, corresponding to boreal 
forests (comparatively larger, broader fine mode particles, with weaker and nearly 
spectrally-neutral absorption) and grass/shrub burning (smaller, narrower particles with 
stronger absorption, becoming more strongly-absorbing as wavelength increases from 
440 nm to 1020nm). Mongu in the southern African savannah exhibits even stronger 


c n 
o 
c 

C/) 

C/) 


TJ 

CD 

"D 

CD 


( / > 
O 
C 

c n 
c f > 


"0 

CD 

T3 

CD 


ACPD 

13, 25013-25065, 2013 

Smoke aerosol 
properties 

A. M. Sayer et al. 


Title Page 


Abstract ■ Introduction 


Conclusions ■ References 


Figures 


C/5 

o 

c 

< 7 ) 
( / ) 


TJ 

CD 

"O 

CD 



Full Screen / Esc 


( / > 
o 
c 

C/3 

GO 

o' 

13 

TJ 

CD 

"O 

CD 


Printer-friendly Version 


Interactive Discussion 


25035 

















absorption, with average SSA around 0.85 in the midvisible. These families can serve 
as candidate sets of aerosol microphysical/optical properties for use in satellite AOD 
retrievals, which are reliant on assumptions about aerosol properties due to the limited 
information content available from existing passive spaceborne imaging radiometers. 
This does not, however, alleviate the difficulty of assuring that an appropriate micro- 
physical model is used for any particular individual pixel-level satellite retrieval. The 
similarity between near-source and transported smoke properties implies that these 
models can be used for AOD retrieval over ocean, as well as land. These models 
represent the climatological properties of biomass burning at these sites, and cannot, 
however, capture the full range of variability at a given site. 

An important outcome of this analysis is that the microphysical models adopted 
by commonly-used satellite AOD retrieval algorithms over ocean are insufficiently ab- 
sorbing to represent aerosol optical properties well in many of these biomass-burning 
regimes. A consequence of this is that they are likely to underestimate the AOD in 
some smoke outflow regions. As these satellite datasets are increasingly used in cli- 
mate applications, and as an evaluation tool for chemistry transport models, this is 
potentially a significant shortcoming. Rare cases of Sahelian smoke suggest that even 
stronger absorption may be seen (Johnson et al., 2008; Eck et al., 2010), which would 
further exacerbate these biases (although in this particular situation the smoke is al- 
most always mixed with more weakly-absorbing dust). Until the launch of future satel- 
lite sensors with increased measurement capabilities, it is important that the continual 
evolution of algorithms using existing sensors includes the adoption of more realistic 
aerosol microphysical models as our knowledge of aerosol properties increases. 
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Table 1. Geographical/sampling information for biomass-burning AERONET land sites, along 
with key references for burning in these regions. 
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Table 2. Size distribution parameters for biomass burning aerosols from the climatological 
AERONET smoke sites (relative to fine/coarse AOD at 550 nm). Note that the fine-mode vol- 
ume relationship is calculated for r f 550 = 0.5 as in Fig. 4 and so will vary for different AOD. r 
indicates Pearson’s linear correlation coefficient. 
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0 -2207,550 
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1 -57 T c,550 

3.34 

0.67 

Mukdahan 

1.51 

0. 1 82t, 550 
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Table 3. As Table 2, except expressed in terms of AOD at 440 nm. 
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0.137T f 440 

0.138 +0.010T f440 , r = 0.20 

0.362 +0.017t i44 o, t = 0.14 

1 .29T C 440 

2.80 

0.69 

Yakutsk 

3.89 

0 . 120 r f 440 

0.164 +0.01 It, 440 , r = 0.32 

0.461 +0.030t i44 „, r = 0.31 

1 .53t c 440 

3.34 

0.72 


o 


C /) 

o 

c 

</) 

C/) 


"0 

■o 

CD 


A. M. Sayer et al. 


Title Page 


■ 


Abstract 


Introduction 




Conclusions 


References 


u 

(/)' 

o 

c 

C/) 

Cfl 

o' 


"0 

0) 

“O 

CD 


D 

c/)‘ 

o 

c 

Cfl 

(J) 

o’ 


TJ 

0) 

"O 

CD 


Figures 



Full Screen / Esc 


Printer-friendly Version 
Interactive Discussion 



25048 

















ACPD 

13, 25013-25065, 2013 


g 

C/5* 

o 
c 
c n 
CO 
o' 


"0 


Q) 

T3 

CD 


Smoke aerosol 
properties 


D 

(/)' 

o 

c 

Table 4. Refractive index and SSA for biomass burning aerosols from the eight near-source 
AERONET sites. Note that SSA is calculated for r f 550 = 0.5, r c550 = 0.03 as in Fig. 4 and will 
vary for different AOD. 
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440 nm 

Refractive index (n-/'k) 
675 nm 870 nm 

1020nm 

440 nm 
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675 nm 870 nm 

1020nm 


Alta Floresta 

1 .46— 0.01 1i 

1 ,48-0.0094i 

1 ,48-0.0085i 

1 ,47-0.0082i 

0.92 
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Table 5. Geographical/sampling information for AERONET coastal/island sites exhibiting cases 
of transported smoke. 


Site 

Latitude, 

degrees 

Longitude, 

degrees 

Elevation 

m.a.s.l. 

Number of 
retrievals 

Dates (YYYYMMDD) 

Ascension Island 

-7.97600 

-14.4150 

30 

19 

20030902-06, 20080904, 20080906-07, 
20080919, 20080925, 20080928 

Bach Long Vy 

20.1330 

107.733 

5 

10 

201 00506, 201 00509, 201 1 0209, 201 1 0224, 
201 1 031 4, 201 1 041 8, 201 1 0421 , 201 1 0423 

Barrow 

71.3120 

-156.665 

0 

5 

20040703, 20100612 

CEILAP Buenos Aires 

-34.5670 

-58.5000 

10 

28 

20010803, 20010921-22, 20040813-15, 
20040827, 20040904, 20060722, 20060824, 
20060910, 20060922, 20100821-22, 20110903 

COVE 

36.9000 

-75.7100 

37 

12 

20020706-09 

Darwin 

-12.4240 

130.892 

29 

20 

20071013, 20071016, 20071018-19, 
20091011-12, 20091104 

Hornsund 

77.0010 

15.5600 

10 

6 

20060502-03 

Inhaca 

-26.0410 

32.9050 

73 

33 

20000823, 20000831 , 20000902-03, 
20000905-06, 20000910, 20000914, 
20001004-05, 20001007, 20001009-10, 
20010820-21, 20010823, 20010829, 
20010908,20010910-11 20010916 

La Jolla 

32.8700 

-117.250 

115 

7 

20031028, 20070804, 20071024-25 

Monterey 

36.5930 

-121.855 

50 

4 

20080627, 20080710, 20080712 

Noto 

37.3340 

137.137 

200 

5 

20030606, 20080422, 

San Nicolas 

33.2570 

-119.487 

133 

2 

20031027 

Saturn Island 

48.7830 

-123.133 

200 

22 

20080630, 20080701-02, 20100802, 
20100804, 20100806, 20100816-17 

Sevastopol 

44.6160 

33.5170 

80 

28 

20070809-10, 20070901-03, 20080801, 
20100815-18 

Singapore 

1 .29800 

103.780 

30 

10 

20090806-07, 20110905-06, 20120924 

Trinidad Head 

41.0540 

-124.151 

105 

13 

20060925-26, 20080709 

UCSB 

34.4150 

-119.845 

33 

19 

20031025-27, 20070817-18, 20071022-23, 
20071025-26 
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Fig. 1. Locations of AERONET sites used in this work. Red diamonds indicate the near-source 
sites, and blue triangles the coastal/island sites with occassional cases of transpored smoke- 
dominated aerosols. 
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Fig. 2. Median retrieved size distributions from biomass burning cases for near-source 
AERONET sites considered (Fig. 1, Table 1). Figures in parentheses indicate the number of 
retrievals for each site/geographical region. 
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Fig. 3. Example of linear fitting procedure used to arrive at relationships between AOD and 
fine-mode microphysical model parameters, for Alta Floresta. Black diamonds indicate points 
used in the fit (green line), while excluded outliers are shown with red triangles, (a) shows fine 
mode volume radius (r vf ), and (b) the fine mode spread (<j f ). 
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Fig. 4. Properties of smoke aerosols for near-source sites, for a reference fine-mode AOD of 0.5 
and coarse-mode AOD of 0.03 at 550 nm. Also shown is the “fine-dominated” model of Sayer 
et al. (2012a) used in the SeaWiFS Ocean Aerosol Retrieval (SOAR) dataset. Panels show 
(a) size distributions, and spectral (b) extinction, (c) single scatter albedo, and (d) asymmetry 
parameter. 
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Fig. 5. Annual cycle of AERONET-retrieved aerosol properties and water vapor at Bonanza 
Creek and Yakutsk. Panels show (a) r f 550 , (b) columnar water vapour, (c, d) fine mode size, 
and (e, f) aerosol SSA. 
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Fig. 6. Linear relationships between fine-mode median radius and spread, (a, b) consider near- 
source land AERONET sites for volume and number radii respectively (legend in a), and (c, d) 
coastal/island transported smoke cases for volume and number radii (legend in c). 
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Fig. 7. True-colour image from MODIS aboard the Aqua satellite showing smoke (grey-brown 
hues) blowing from fires in California (USA) into the Pacific Ocean on 27 October 2003. Red di- 
amonds show, from North to South, the locations of the UCSB (34.4° N, 1 1 9.8° W), San Nicolas 
(33.3° N, 1 1 9.5° W), and La Jolla (32.9° N, 1 1 7.3° W) AERONET sites. 
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Fig. 8. True-colour image from MODIS aboard the Terra satellite showing smoke (grey-brown 
hues) spreading across South America and into the southern Atlantic Ocean on 22 August 
2010. Image acquired from two consecutive Terra orbits. Red diamonds show, from North to 
South, the locations of the Alta Floresta (9.87° S, 56.1°W), Ji-Parana SE (10.9° S, 61. 9° W), 
CUIABA-MIRANDA (1 5.7° S, 56.1° W), and CEILAP-Buenos Aires (34.6° S, 58.5° W) AERONET 
sites. The green line shows the HYSPLIT 10-day back-trajectory for the air mass ending at 
2 km above Buenos Aires at 00:00 UTC on 23 August 2010; triangles indicate the position at 
00:00 UTC each day. 
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Fig. 9. Time series of aerosol properties at AERONET sites shown in Fig. 8 for 10-30 August 
2010. Panels show (a) t 440 , (b) a (both from the direct-Sun AERONET observations), (c) r eff f , 
and (d) SSA at 440 nm (both from the AERONET inversion dataset). The shaded grey area 
indicates 21-23 August when the smoke plume was observed at Buenos Aires. 
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Fig. 10. Properties of smoke-dominated aerosol from near-source sites located in the Amer- 
icas (Fig. 2), and transported smoke cases for which these sites provide the closest match. 
From left-right, plots compare fine mode radius, fine mode spread, aerosol SSA (different sites 
offset slightly for clarity), and coarse mode properties. Grey triangles show the data (with out- 
liers removed, as described in the text) and solid lines the linear fit (Fig. 3). In the SSA plots, 
the shaded grey area indicates the central 68% of retrieved values at this site. Dashed lines 
indicate properties for unpolluted marine aerosol from Sayer et al. (2012b); in the coarse-mode 
properties plot, the dashed box encloses the range of average radii/spreads for marine aerosol 
for different wind speed regimes. In all plots, coloured symbols indicate the cases of transported 
smoke from different sites, as indicated by the legend. 
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Fig. 11 . As Fig. 10, except for near-source sites located in Africa (Fig. 2). 
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Fig. 12. As Fig. 10, except for near-source sites located in Asia (Fig. 2). 


25062 


Coarse spread (o c ) Coarse spread (o c ) 


ACPD 

13, 25013-25065, 2013 


1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 

1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Coarse volume radius (r vc ), pm 

1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 

1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Coarse volume radius (r v>c ), pm 


(h) Yakutsk 



o 


(d) Mukdahan 



O 


( f > 
O 

c 
c n 
CO 
o' 


"0 


Q) 

"O 

CD 


D 


c n 
o 
c 
<y) 
c n 


TJ 

O 

"O 

CD 


Smoke aerosol 
properties 

A. M. Sayer et al. 


Title Page 


■ 


Abstract 


Introduction 




Conclusions 


References 


g 

c ri 
o 
c 
V) 
V) 

o' 

o 

"0 

"O 

CD 


D 

( / )’ 
o 
c 
c n 

(j) 

o' 

13 

"0 

"O 

CD 


Tables H Figures 



Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 

















uiri ‘( J A -0 snipej omnjoA ouij 


ACPD 

13, 25013-25065, 2013 


(a) Jabiru 



550 nm fine AOD 


(b) Jabiru 



550 nm fine AOD 


(c) Jabiru 



l.O 
0.9 

2 0.8 

1 

§4 0.7 

1 0.6 
U 

0.5 
0.4 

1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Coarse volume radius (r vc ), Jim 


(d) Jabiru 



Fig. 13. As Fig. 10, except for near-source sites located in Oceania (Fig. 2). 
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Fig. 14. Error on retrieved AOD at 550 nm as a function of AERONET AOD at 550 nm, from 
MODIS, SeaWiFS, and MISR data products at Ascension Island. Coloured lines indicate the 
least-squares linear fit of bias vs. AOD for each dataset. The number of matches for each 
sensor is shown in parantheses. 
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Fig. 15. TOA reflectance and AOD retrieval errors for four wavelengths, (a-d) show TOA re- 
flectance as a function of t 550 , for strongly (black), moderately (red), and weakly (blue) ab- 
sorbing aerosol microphysical models over ocean, (e-f) show the error in AOD which would be 
retrieved if measurements at that wavelength were used and a moderately- or weakly-absorbing 
aerosol assumed, if the true aerosol were instead strongly-absorbing. 
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